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THE SOLAR SIBLINGS IN THE GAIA ERA 

C.A Martmez-Barbosa^, A.G.A. Brown^ and S. Portegies Zwart^ 


Abstract. We perform realistic simulations of the Sun’s birth cluster in 
order to predict the current distribution of solar siblings in the Galaxy. 
We study the possibility of finding the solar siblings in the Gaia cat¬ 
alogue by using only positional and kinematic information. We find 
that the number of solar siblings predicted to be observed by Gaia will 
be around 100 in the most optimistic case, and that a phase space only 
search in the Gaia catalogue will be extremely difficult. It is there¬ 
fore mandatory to combine the chemical tagging technique with phase 
space selection criteria in order to have any hope of finding the solar 
siblings. 


1 Introduction 

It is commonly accepted that disk stars are born in open clusters (Lada et al. 
II993|I . Particular features of the solar system suggest the Sun was born in a 
cluster. For instance, the radioactive elements found in the meteorite fossil record 
(Looney et al. 120061) can be explained by a supernova explosion in the vicinity of 
the newborn Sun. The coplanarity of the Solar system planets and the eccentric 
orbits of the Kuiper-belt objects suggest that the Sun had a close encounter with 
another star belonging to the same open cluster (Morbidelli & Levinson I2004p . 
Such an encounter is expected to have occurred in the early history of the solar 
system IMalhotra I2008|l . 

If the Sun was born in an open cluster 4.6 Gyr ago, the solar siblings -stars 
born in the same cluster- would be spread out over the Galactic disk. The iden¬ 
tification of at least a small fraction of these stars would enable us to understand 
the conditions in which the solar system was formed. 

In this paper we study the evolution and disruption of the Sun’s birth cluster 
in the Galaxy by means of state-of-the-art simulations, taking into account the 
effects of the spiral arms and bar of the Galaxy as well as the internal processes 
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in the cluster. We aim to predict the current phase-space distribution of the Solar 
siblings and the probability of detecting them in the large Gaia database. 

2 Simulations set up 

The Milky Way is modelled as an analytical potential composed by an axisym- 
metric component together with a rotating central bar and spiral arms. We use 
the potential of Allen & Santillan (I1991D to model the axisymmetric component 
of the Galaxy. The central bar is modelled with a Ferrers potential (Ferrers I1877|) 
and the spiral arms are modeled as perturbations of the axisymmetric potential 
(Antoja et al. I201ip . The bar and spiral arms are non transient structures that 
rotate as rigid bodies with different pattern speeds. The Galactic parameters used 
in the simulations are taken from Martmez-Barbosa et al. (|2015ll . 

The Sun’s birth cluster on the other hand is modelled as a spherical distribution 
of stars that obey a Plummer potential (Plummer I191ip . We use a Kroupa initial 
mass function (IMF) (Kroupa I200ip to model the mass distribution of the Sun’s 
birth cluster. The minimum and maximum stellar masses are 0.08 Mq and 100 Mq 
respectively. Additionally, the metallicity of the Sun’s birth cluster is assumed to 
he Z = 0.02. Following Portegies Zwart ((200911 . we set the initial mass of the Sun’s 
birth cluster between 500 and 3000 Mq. The initial radius is set between 0.5 and 

3 parsec. 

The initial phase-space coordinates of the centre of mass of the Sun’s birth 
cluster (xcmWcm) are computed by integrating the orbit of the Sun backwards in 
time, taking into account the uncertainty in its current Galactocentric position 
and velocity. By using this technique, Martmez-Barbosa et al. (I2015II found the 
distribution of possible positions and velocities of the Sun at its birth (P(i?b), 
P{Vh)). Given that for most of the Galactic parameters, (P(i?b)) is peaked around 
9 kpc, we assume initially for the Sun’s birth cluster Ccm = 9 kpc. The Sun’s birth 
cluster is then evolved under the influence of its self-gravity and the strong external 
tidal field produced by the Galaxy. We also consider stellar evolution effects. We 
used the huayno code (Pelupessy et aL I2012p to compute the internal gravitational 
force in the Sun’s birth cluster. To compute the force due to the external Galactic 
potential, we use a high-order Rotating bridge (Martmez-Barbosa et al. 120151) . 
The stellar evolution effects were modelled with the SEBA code (Portegies Zwart & 
Verbunt 119961 Toonen et a,l. I2012|) . These codes were coupled through the AMUSE 
framework (Portegies Zwart et al. 120131) . 

3 Results 

Since the age of the Sun is approximately 4.6 Gyr (Bonanno et aL I2002D . we evolved 
the Sun’s birth cluster over that period of time in the Galactic model described 
previously. We used different initial masses and radii and different combinations 
of bar and spiral arms parameters. At the end of each simulation, the Sun’s birth 
cluster is completely disrupted, with the solar siblings located in different regions 
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Fig. 1. Left: Parallax Vs. proper motion of the solar siblings. Right: distribution of 
radial velocity of the solar siblings. Top: Astrometric properties when the Galaxy has 
two spiral arms. Bottom: Astrometric properties when the Galaxy has four spiral arms. 
We varied the mass and pattern speed of the bar and the amplitude and pattern speed 
of the spiral arms. 


of the Galactic disk. The radial and azimuthal dispersion of the hnal (present day) 
distribution of solar siblings depends on the combination of Galactic parameters. 
From the galactocentric distribution of solar siblings, we compute their parallaxes, 
proper motions and radial velocities. An example is shown in Fig. [l] Here we 
have superposed the astrometric properties of the solar siblings for 130 different 
combinations of bar and spiral arms parameters. The initial mass and radius 
of the Sun’s birth cluster were fixed to 500 Mq and 0.5 parsec respectively. As 
can be observed, the solar siblings span a large range of parallaxes and proper 
motions (e.g. left panel Fig. [^. The dark blue region corresponds to the area 
in phase-space where most of the solar siblings are located. The parallax of the 
solar siblings in this region is in the range 0.3 mas; therefore, most of 

the solar siblings are relatively far from the Sun. Note that the proper motion 
of the solar siblings there is around 6 mas yr“^, which can be explained as stars 
moving along the solar circle with velocities comparable to that of the Sun (Brown 
et al. I2010|) . The distribution of radial velocities is centered around 0 kms“^ (e.g. 
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right panel of Fig. [^. The small relative motion of the solar siblings with respect 
to the Sun is due to the Galactic model used in the simulations. Non-transient 
spiral structure does not produce significant radial migration. Transient spiral 
arms and the interaction of the Sun’s birth cluster with giant molecular clouds 
(GMCs) would produce a broader range of phase space coordinates and therefore 
broader radial velocity distributions. 

Even though in the simulations the Sun’s birth cluster is not much affected by 
radial migration, we found that the expected number of solar siblings predicted 
to be observed by Gaia is not higher than 100. If we assume four spiral arms in 
the Galactic model, this number is reduced to 50 stars on average. Given that the 
Gaia catalogue will contain one billion stars, it will be extremely difficult to find 
solar siblings using only phase space, even if we look for a contrast between the 
simulated phase-space of solar siblings and the disk stars of the Gaia data base (as 
done in Brown et al. (|2010|) 1. Therefore, in order to find possible solar siblings, 
it is necessary to make use of the chemical tagging technique combined with the 
phase space selection criteria established through the simulations. 
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